Manipulating 4/ quadrupolar pair interactions in TDB2C2 using a magnetic field 
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Resonant soft x-ray Bragg diffraction at the Tb Af^s edges and non resonant Bragg diffraction 
have been used to investigate orbitals in TbB2C2- The Tb 4/ quadrupole moments are ordered 
in zero field below Tjv and show a ferroquadrupolar alignment dictated by the antiferromagnetic 
order. With increasing applied field along [ffO] the Tb 4/ magnetic dipole moments rotate in a 
gradual manner toward the field. The quadrupole moment is rigidly coupled to the magnetic moment 
and follows this field-induced rotation. The quadrupolar pair interaction is found to depend on the 
specific orientation of the orbitals as predicted theoretically and can be manipulated with an applied 
magnetic field. 

PACS numbers: 71.70.Ch; 75.40.Cx; 78.70.Ck 



Correlation between conduction electrons and elec- 
tronic orbitals leads to interesting materials properties 
such as metal-insulator transitions, colossal magneto re- 
sistance and superconductivity. Aspheric electronic or- 
bitals, characterized by their quadrupole moment, may 
order and cause partial charge localization of the con- 
duction electrons^ or mediate coupling between cooper 
pairs. 2 Orbital order in / electron materials is domi- 
nated by coupling with the lattice (Jahn- Teller) or by 
indirect Coulomb interactions via the conduction elec- 
trons. That the latter can be important for intermetallic 
compounds was established theoretically 3 ^ but detailed 
experimental knowledge is limited due to the difficulty 
observing the associated orbital excitations. A recent 
neutron diffraction study provided evidence of a modu- 
lated quadrupolar motif in PrPb3, believed to be a direct 
consequence of indirect Coulomb interactions exhibiting 
oscillatory nature^ The large orbital momentum of the 
/ electronic shell also gives rise to a significant influence 
of higher multipole moments and may lead to hidden or- 
der phase transitions as demonstrated in the extensively 
studied URu2Si2^ Therefore it is important to under- 
stand the quadrupolar and higher order multipole pair 
interactions in these materials. 

Quadrupolar order has been successfully investigated 
using neutron scattering in applied fields where induced 
magnetic moments reveal the underlying quadrupolar ar- 
rangement or motif. In addition, the relatively weak x- 
ray diffraction intensity of the quadrupole moments can 
be observed using synchrotron techniques allowing a di- 
rect determination of orbital motifs. Resonant x-ray scat- 
tering at the L2,3 edge provided the first proof of the or- 
bital motif in DyB2C2 TbB2C2 is proposed to exhibit 
a transition from antiferromagnetic (AFM) to antiferro- 
quadrupolar (AFQ) order in an applied magnetic field 9 
and is therefore an interesting candidate for investigation 
of its orbital interactions. 

A magnetic field is time-odd and cannot couple to 
quadrupole moment which is time-even, but neverthe- 
less in TbB2C2 the ordering temperature increases to 



35 K, well above Tjv = 21.7 K, in an applied magnetic 
field of 10 T. Consequently the interplay between dipole 
(time-odd) and quadrupole (time-even) moments can be 
readily investigated. Neutron diffraction studies revealed 
similarities between the magnetic structure in applied 
fields and the magnetic structure in the combined AFQ 
and AFM phase of DyB2C2. 10 However, resonant x-ray 
diffraction at the Tb L 3 edge suggests AFQ is present be- 
low Tjv in zero field with a k vector component of (0,0, i) 
as found for DyB2C2iii 

To test the suggested field-induced quadrupolar order- 
ing we have investigated TbB2C2 in an applied mag- 
netic field, H, with x-ray scattering techniques. We 
show that orbital order is present in the AFM phase 
of TbB2C2 at zero field, however, with ferroquadrupolar 
alignment which is dictated by the AFM structure. The 
quadrupoles are rigidly linked to the magnetic dipole mo- 
ments and rotate with H along the [110] direction. We 
show that, due to this rotation, the quadrupole pair in- 
teractions become stronger with applied field. 

A single crystal of TbB2C2 was grown by Czochralski 
method using an arc-furnace with four electrodes. Sam- 
ples were cut and polished perpendicular to the [001] di- 
rection. The (00^) reflection was recorded at the Tb 
M 4i5 edges of TbB 2 C 2 at the RESOXS end-station of 
the SIM beamline at the Swiss Light Source of the Paul 
Scherrer Institut. A permanent magnet provided a field 
of 1 T parallel to the [110] direction in the scattering 
plane. In addition, non resonant x-ray Bragg diffraction 
of the (01 1 ) and (11 |) reflections was performed at the 
BL19LXU beam line at SPring-8 with 30 keV x-rays us- 
ing a Ge solid-state detector to eliminate higher order 
harmonics. A cryomagnet provided H along the [110] 
direction. In both sets of measurements, structure fac- 
tors were derived from the integrated intensity, corrected 
for polarization, Lorenz factor, absorption and sample 
geometry. 

The resonant Bragg diffraction results, Figure^, show 
the energy dependence of the (00^) reflection recorded 
at the Tb edges in zero field. The energy profile 
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FIG. 1: Energy profile of the (00±) reflection at the Tb M 4 , 5 
edges in TDB2C2 in (a) zero field and (b) applied field of 1 T 
along [110], recorded with a or n incident radiation at 11 K. 

is independent of sample rotation about the Bragg wave 
vector. A significant part of the diffracted intensity is 
due to the total reflectivity of the polished sample. The 
(00 i) diffracted intensity is relatively weak and appears 
magnetic in origin as witnessed by the relatively large 
intensity for incoming x-rays polarized in the scatter- 
ing plane (tt). In contrast, the diffracted intensity with 
H = 1 T, along the [110] direction, shown in Fig. ^ is 
much stronger and strongest for incoming radiation po- 
larized perpendicular to the scattering plane (er). The 
energy profiles shown in Fig. ^ are expected to be in- 
sensitive to the magnetic structure, not to be confused 
with the total intensity that does depend on the mag- 
netic structure. The drastic change with H indicates the 
presence of either an additional component of resonant 
scattering or a change in relative contribution of two dif- 
ferent components of resonant scattering. The energy 
profile recorded in an applied magnetic field is similar to 
that recorded for the (00^) reflection in DyB2C2 in the 
AFQ phase. 12 This implies that for TbB 2 C2, under the 
influence of an applied magnetic field of 1 T, the (00^) 
reflection is dominated by quadrupolar scattering. 

The non resonant Bragg diffraction results illustrate 
that the quadrupole moment follows the magnetic mo- 
ment rotation in applied fields. Fig.[2K,b shows the struc- 
ture factors that relate directly to the order parameter, 
obtained from the integrated, non-resonant Bragg inten- 
sities of the (01 1 ) and (ll^r) reflection as a function 
of applied field. In the AFM phase, at 15 K, the struc- 
ture factors show a linear increase for applied fields below 
1.5 T. In the paramagnetic phase, at 27 K, the ordered 
phase is entered above 3 T. The (0l|) reflection appears 
at 3.2 T while the (H^) reflection appears at 3.5 T. 
However the evolution of the structure factor with in- 
creasing field is smooth for both reflections without a 
hint of a phase transition at 3.5 T for (0l|). A sim- 
ilar behavior is observed as function of temperature in 
an applied field of 5 T (not shown). The intensities of 
(0l|) and show a smooth decrease with increas- 
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FIG. 2: Structure factor of the (01§) and (11^) reflection as 
function of applied magnetic field recorded with 30 keV at (a) 
T = 27 K and (b) T = 15 K. (c) and (d) show the deduced 
angles of Tb 4/ orbital rotation as defined in Fig. [3] 

ing temperature and disappear below 30 K and 29.5 K, 
respectively. 

A structural transition with alternating displacements 
of the B and C atoms gives intensity in both type of 
reflections^ and cannot account for the present observa- 
tion of intensity at (0l|) and the absence of intensity at 
(ll4r). Magnetic, or time-odd, x-ray diffraction is gen- 
erally weak and does not explain the magnitude of the 
observed intensity. We conclude that the scattering arises 
from time-even x-ray diffraction, i.e. from the aspheric 
charge distribution of the 4/ shell, as such scattering can 
be intense. 

To determine the structure factor it is assumed that 
the Tb ion on each site has the same charge distribution 
except for its specific orientation. Using the formalism for 
non resonant time-even scattering of Lovesey et. alm^ we 
deduce two independent Tb ions contribute and obtain 
for the structure factors: 

F n i cx -2 [sin(20i) - sin(20 3 )] (Q u « ), (1) 
F 0l i cx [008(20!) - cos(20 3 )] (Q 01 i ), (2) 

when the Miller index I equals an odd integer. The prin- 
cipal axis of the 4/ orbital at atom n is canted in the 
afe-plane with respect to the [100] direction with angle 
4>n- (Qhki) is the expectation value of the Tb 4/ time- 
even multipolar moment for a given field, temperature 
and Bragg wave vector k = (k a ,ki,,k c ): 
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(Q ui ) cx ^ [(j 2 )*j + 3V3<j 4 >*2 - \/l82(j 6 )*!j , (3) 

(Qoii)«<Q lli ) + ^^l(j 4 )*t. (4) 

where is a structure factor of the chemical unit cell 
that is a linear sum of the atomic tensors, ^=(T^) — 
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FIG. 3: (a) Phase diagram of TbB2C2 with magnetic field applied along [110]. (b) and (c) illustrate the orientation of magnetic 
moments (arrows) and quadrupole moments (ovals) in TDB2C2 in phase Ha in zero field and under applied fields, respectively, 
(d) and (e) illustrate phase lib under the same conditions. 



(T* q ). describes the asphericity of the Tb 4/ elec- 

tronic shell where (Tf) equals the quadrupole moment, 
(T2) the hexadecapole moment and (Tf) the hexacon- 
tatetrapole moment. (72), (J 4} and (je) are Bessel func- 
tion transforms of the radial distribution of the 4/ elec- 
tronic charge. In eqs.03and21 (k c /k) 2 ~ 1 and we neglect 
higher order terms (kb/k) 4 and (kb/k) 6 as kb/k <Cl. 

We discuss the configuration of the quadrupoles but 
the order of the higher multipoles may alter also. When 
the 4/ quadrupoles are aligned along [110] (</>i=</>3=|7r) 
as drawn in Fig. |3fb), F Q1 ^ and F n ^ are both zero. 
When the 4/ quadrupoles tilt toward the applied field 
with angle S (<f>i=jn + S; 3 =f7r-<5), Fig.|3Jc), F 01 i is 
finite and F-, , ; remains zero. Neutron diffraction shows 
that in the AFM phase the dipole moments tilt away from 
(110) with angle A, Fig.^d) (^i=f tt-A; (p3=jir-A)M 
In this case, both F 01 i and F xl i are zero in zero field 
while under an applied field both become finite (Fig.^e), 
&=f 7T - A + 6; 03=|7r - A - 6). 

Consequently, at 27 K, between H = 3.2 and 
H = 3.5 T, the magnetic and quadrupolar structure of 
TbB2C2 is characterized by phase Ha. The quadrupoles 
are aligned parallel to [110] as illustrated in Fig.|3{b) and 
are tilted in applied field as illustrated in Fig.|3fc). Above 
3.5 T the compound is characterized by phase lib. The 
quadrupoles are tilted away from the [110] direction and 
neighboring atoms in the afr-plane are tilted in opposite 
directions as illustrated for H = in Fig. |3Jd) and for 
H 7^ in Fig.|sje). We show the phase diagram schemat- 
ically in Fig.^a) and define a reorientation temperature 
Tr by the boundary between phases Ha and lib. 

Multipole moment rotation in an applied field is linear 
for small 5, A (see Eq. Q) hence the linear increase of 
the structure factor as observed at 15 K (see Fig. Efb)). 
The magnetic moment of the 4/ shell tilts toward the 
direction of the magnetic field and the time-even mul- 
tipole moments follow accordingly, giving the observed 
contrast. This shows that (i) the coupling between mag- 
netic and quadrupole moment is rigid and (ii) the orbitals 



are already ordered in zero field, although with a ferro- 
quadrupolar alignment. This is consistent with the result 
in Fig.^, as the ferroquadrupolar alignment in zero field 
does not contribute to scattering at (00^), whereas in 
an applied field, due to the orbital rotation, time-even 
scattering dominates (Fig. [IJa) . 

The normalized structure factor of (0l|) and (11^) 
reflections at 5 T with 9 < I < 17 and are presented in 
Fig 21 They are equal within the experimental uncer- 
tainty which shows that ^\ is small (Eqs.lSjand^J. The 
solid line corresponds to Eqs. Q and El with <E , 2/'^ , ! = Q-5 
and ^2/^1=0. 5 and shows that the intensity can be at- 
tributed to time-even x-ray diffraction. The limited range 
of k makes the determination of ty\ uncertain and this 
analysis merely confirms the absence of scattering from 
B and C atoms in contrast to DyB2C2>^ Thus orbital 
order is present independent of a structural transition 
in TbB 2 C2 verifying the orbital interaction is mediated 
via the conduction electrons. A similar conclusion was 
reached for DyB2C2 from inelastic neutron scattering^ 

With (Q 01 i) = (Qui) for a given field and tempera- 
ture the angle A is deduced from the ratio between F 01 a 
and -fii^s . In addition the quantity sm(2d) (Qhki) follows 
from Eqs. n an d El and S has been estimated assuming 
(Qhki) is constant as a function of field and temperature. 
Both results are presented in Figs.|2tc) and (d). If (Qhkl) 
increases with H the actual values for S are lower than 
presented in Fig. Nevertheless, the increase in 5 is 
gradual as a function of the applied field and reflects the 
orbital moment rotation. 

The angle A is a result of the competition between 
the indirect Coulomb interaction and the magnetic ex- 
change interaction between neighboring ions in the ab- 
plane. The first interaction favors perpendicular align- 
ment of the quadrupoles while the latter favors a parallel 
alignment. An increase in A is observed with H. This 
shows that the relative strength of the indirect Coulomb 
interaction increases as the quadrupoles rotate perpen- 
dicular to each other in an applied field. To the authors 
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FIG. 4: Structure factor of (01~) and (lli) reflections at 
B=5 T and T=2 K. The solid line is a fit to Eq.0 

knowledge, this is the first time that an increase in indi- 
rect Coulomb interaction has been observed in a straight- 
forward manner. We also note that the determination of 
A is robust against changes in (Qhki) with applied field. 

The increase in ordering temperature with applied 
fields is in line with an antiferroquadrupolar interaction 
which becomes stronger when the angle between the or- 
bitals moves from parallel to perpendicular alignment. 
The magnetic moments interact via the polarized spins of 
the conduction electrons while the multipolar moments 
interact via the charge density of the conduction elec- 
trons. In other words, charge screening by the conduc- 
tion electrons changes when the 4/ orbitals rotate and 
consequently the magnitude of the quadrupole pair inter- 
actions changes. This is consistent with the theoretical 
framework of Teitelbaum and Levy^ where electric mul- 
tipole coupling via the conduction electrons depends on 
the relative angular positions of the ion cores. 

The occurrence of combined AFQ and AFM below T/v 
is consistent with the quasi doublet ground-stated Re- 



moval of the degeneracy of a doublet ground-state may 
account for one phase transition. Relaxation between 
the two singlets takes place above T/v and the magnetic 
and quadrupolar moments fluctuate. The increased en- 
ergy separation between the two singlets below T/v is due 
to the magnetic exchange interaction, but both the mag- 
netic and quadrupolar moment of the ground state singlet 
are observed, and this is illustrated by this study. 

In conclusion, no magnetic field induced ordering of 
quadrupoles exists in TbB2C2- The orbitals are already 
ordered in the AFM phase in zero field with a ferro- 
quadrupolar motif. On increasing the applied field along 
[110] the Tb 4/ magnetic moments rotate toward the field 
direction in a gradual manner. The quadrupole moment 
is rigidly coupled to the magnetic moment and follows 
this rotation. Neighboring Tb 4/ orbitals move from 
parallel to perpendicular alignment and the quadrupo- 
lar pair interaction increases, witnessed by an increase of 
angle A between neighboring quadrupoles. This study 
shows that the quadrupolar pair interaction depends on 
the specific orientation of the orbitals as predicted for in- 
direct Coulomb interaction via the conduction electrons 
and can be manipulated with an applied field. 
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